Generic features of the spectrum of trapped polarized fermions 
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We show that bimodal radio frequency spectra universally arise at intermediate temperatures in 
models of strongly interacting trapped Fermi gases. The bimodality is independent of superfluidity 
or pseudogap physics, depending only on the functional form of the equation of state - which is 
constrained by dimensional analysis at low temperatures and the virial expansion at high tempera- 
tures. In addition to these model independent results, we present a simple calculation of the radio 
frequency line-shape of a highly polarized Fermi gas which uses energetic considerations to include 
final state interactions. While this model only qualitatively captures the line-shapes observed in 
the experiments, it provides a conceptually clean and powerful technique for estimating the energy 
scales and how they vary with experimental parameters. 
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Can spectroscopy be used to detect pairing in a gas 
of fermionic atoms? The paradigm for thinking about 
this question was set in 2003, when Greiner, Regal and 
Jin [l| presented an experimental study of the microwave 
spectrum of a two component gas of 40 K on the BEC 
side of resonance (where two potassium atoms are capa- 
ble of forming weakly bound molecules). Their sample 
consisted of a noncondensed gas of diatomic molecules in 
chemical equilibrium with a gas of atoms. They found an 
easily interpreted bimodal spectrum: a sharp line came 
from the excitations of free atoms, and a broad peak 
from the dissociation of molecules. Later experiments on 
colder samples in different parameter ranges showed simi- 
lar bimodality and were interpreted in similar ways d El ■ 
In particular, when such a bimodal peak was seen by 
Chin, Bartensein, Altmeyer, Riedl, Jochim, Denschlag, 
and Grimm [2| on the BCS side of resonance, it was 
taken as a sign of Cooper pairing. This interpretation 
was further reinforced by mean field calculations which 
showed that the finite temperature trapped paired gas 
does indeed produce bimodal spectra 0, H, H Q ■ This 
paradigm has been shattered by the realization that there 
exist models which produce bimodal spectra in the ab- 
sence of pairing Here we give a simple argument for 
why trapped strongly interacting Fermi gases generically 
show a bimodal spectrum, irrespective of the presence of 
pairing. 

Our argument relies on two ingredients: a harmonic 
trap plus qualitative features of the equation of state. We 
establish that the equation of state of the unitary Fermi 
gas has these features through a dimensional analysis ar- 
gument supplemented by the first terms of the virial ex- 
pansion. Consistent with previous results 0, d, HUH, 
we conclude that one of the spectral peaks comes from 
atoms at the edge of the cloud, and the other from atoms 
at the center. This is similar to the mechanism by which 
multiple spectral peaks appear in the spectra of clouds 
of Bosons in an optical lattice [1 El • 

We quantify our arguments by presenting a simplified 



calculation of the spectrum of a trapped two-component 
Fermi gas in the limit n \ /nt — > 0. Although we include 
final-state interactions [ll|, El El El El El, El EI, 
Tij l in our quantitative calculations, we emphasize that 
the bimodality is not a final-state effect. In the highly 
polarized limit, final state effects set the frequency scale 
for the spectral line but do not significantly change its 
shape or temperature dependence. Moving away from 
this limit of large polarizations, final state interactions 
can have more profound effects [3, El El El El, El, 

II El, El H El III. 



We consider an experiment on a gas of atoms with 
two spin state: j and J,. We imagine that a probe (for 
example radio waves) excites atoms from the J, state 
to a third state x. In most recent experiments on 6 Li 
the states involved are the three lowest energy hyperfine 
states |T) = |1>, ||) = |2) and \x) = |3), (E 1 < E 2 < E 3 ) 
though other combinations are possible [2^|. The spec- 
trum I{v) measures, for a fixed probe intensity, the rate 
of population transfer from J, to function of the 

detuning v from the free-space resonance. 

In principle, moments of this spectrum can 
be calculated from the sum rule [H El, El, 
v=S<LvvI{v)l Sdvl{v) = ±fd 3 r [v u (r)-v n (r)) 

(ip^(r)ip^(0)ipi(0)'4''\(r)}, where v.ij(r) is the interaction 
potential between atoms in states i and j, and i/v is an 
annihilation operator. As pointed out by Pethick and 
Stoof [231 ] . extracting useful information from this result 
is problematic because this expectation value (which is 
closely related to the expectation value of the interaction 
energy) is not a low energy observable: different poten- 
tials which give rise to the same low energy scattering 
properties will give completely different values for P; 
hard spheres have v — 0, while point interaction yield 
v = oo. These potential dependent features come from 
the largely unobservable ultraviolet tail of the spectrum. 
The peak of the spectral line occurs at relatively small v 
and only depends on low energy properties such as the 
scattering lengths ay. 
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Here, we explicitly restrict our discussion to the peak of 
this extended spectrum, beginning with the conceptually 
simple limit where there is a single down-spin particle 
in a gas of up-spin particles with uniform density nj. 
Barring pathologies in the matrix elements of the probe, 
and subtleties associated with vertex corrections |23| , the 
peak in I(i>) should be located at detuning 



(1) 



where e (T (n-f) is the energy of a single particle of spin a in 
the sea of | atoms. Note that e a is neither the self-energy, 
nor the chemical potential, though it can be derived from 
the self- e nerg y by integrating with respect to the coupling 



constant 



24j. The fact that e CT ^ \i a is clear from con- 
sidering the dilute high temperature limit, where \i a is 
very negative, but e a is very small. Assuming that the 
inhomogeneous broadening dominates the lineshape, one 
can approximate the homogeneous lineshape as a delta 
function. In a trapped gas, applying this approximation 
locally will yield an inhomogeneous spectrum 
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Note that we are explicitly neglecting the role of "Frank- 
Condon" factors, related to overlaps between the states 
in the two channels, and do not include "excitonic" effects 
related to pairing in the final state [H, ES E3, El El • To 
correctly model the uniform spectrum, as experimentally 
studied in 21, 22|, it would be essential to include such 
factors. 

Given that tt.^ is vanishingly small, and that for short- 
range interactions atoms in the | state will not interact 
with one-another, we take n-\{r) to be given by the stan- 
dard non-interacting result. 
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where we have used the Thomas-Fermi approximation 
to account for the slowly varying trap potential V[r). 
For harmonic traps, this Thomas-Fermi approximation 
is excellent, even at zero temperature and for small num- 
bers of particles 25|. The density is expressed in terms 
of the polylogarithm function 33/2(2) = J^. z J /j 3 / 2 and 
the thermal wavelength = 2irh 2 /rnksT. Similarly the 
down-spin density is 
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where we have assumed /i| < and |/3/ijJ 3> 1. Defining 
the function fj,(v) to be the inverse of u p , i.e. Vp\ji(y)] = 
v, and assuming a harmonic trap: V(r) oc r 2 , the spec- 
trum is 



I(y) oc -y/^ 
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FIG. 1: Energy of a single minority species atom in a Fermi 
sea of majority species calculated from Nozieres and Schmidt 
Rink free energy. Inset shows energy scaled by majority 
species density rif and the factor 7 = (ktT) 1 ^ 2 (m/fi 2 ) 3//2 - 
horizontal lines correspond to leading order virial expansion. 
From bottom to top a^/mfcsT/ft 2 = 00, —10, —5, —2. 



We now argue that for short range interactions, e CT (/L<-|-) 
must have the generic form illustrated in figure [TJ By 
considering a virial expansion, we know that for 8a^ suf- 
ficiently negative, e must be proportional to n-\ [261 ] . and 
hence will be exponentially small in /3/if . For sufficiently 
large /3/u-p, the scattering will be unitarity limited and e 
will be linear in /3/i-f. Different theories can only change 
the slope of the latter relationship, and change the ex- 
act location of the crossover between these two limiting 
behaviors. Regardless of these two details, the result- 
ing inhomogeneously broadened spectrum will have the 
following features: (i) At all temperatures, the edge of 
the cloud will contribute spectral weight at v = 0: for 
a harmonic trap a weak \J — \o%{J5v) divergence will be 
found, (ii) At high temperatures < 0, |/3^tl ^ 1)) 

all spectral weight will be in this peak. Including addi- 
tional sources of broadening will obscure the singularity, 
yielding simply an asymmetric peak, (iii) At very low 
temperatures (j3fi 3> 1) only an exponentially small frac- 
tion of the spectral weight will be in this low detuning 
peak, and it will not be observable. Rather, most of the 
spectral weight will come from the center of the cloud 
where most of the down-spin particles reside. This will 
give a single broad peak. At sufficiently high central 
density (n^a 3 3> 1), the location of this peak will scale 
linearly with the central up-spin chemical potential. The 
constant of proportionality will depend on final-state in- 
teractions, (iv) At intermediate temperatures a bimodal 
structure will be seen. Illustrative spectra are shown 
in figure O It should be clear that our interpretation 
of the spectra is vastly different from the canonical one 
0, 0] , where the two peaks are attributed to "atoms" and 
"pairs" . 

Although the qualitative spectral features are generic, 
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FIG. 2: Radio frequency spectra of a harmonically trapped 
gas in the limit of vanishing Ni: number of atoms re- 
maining in the 2 state as a function of the detuning v of 
the probe from the vacuum transition value (offset verti- 
cally for clarity). In all cases a\2 = oo. Top figure mea- 
sures frequencies in physical units for experimentally rel- 
evant parameters: from top to bottom «i3 sj mfcsT/ft 2 = 
-3,-2.6,-2.4,-2.4,-1.8 with T/T F = 1.9,1,0.9,0.7,0.5 
and E f = k B T F = 260/ikHz, 360/ikHz, 360/ikHz, 360ftkHz, 
where Ef = (SNujxUiyUJz) 1 ^ 3 /tt measures the number of ma- 
jority species atoms, with ujj the small oscillation frequency 
in direction j, and h is Planck's constant. Lower figure 
measures frequencies in units of Ef. From top to bottom 
T/Tf = 1.25, 0.7, 0.5, while from left to right a^mk B T/h 2 = 
— 5,-2.55, —1. The shape of the spectrum is mainly set by 
T/Tf, while the frequency scale depends on the interaction 
strength. The upper figures are convolved with a Gaussian 
of width 0.003-E/, while no broadening is used in the lower 
figures. 



we make quantitative predictions by calculating the en- 
ergies in figure [T] from the generalization of the Nozieres 
and Schmitt-Rink free energy [27| to unequal densities. 
Within this approximation, which can be thought of as 
pairwise summing all scattering events, taking into ac- 
count medium effects by a phase space reduction factor, 
the interaction contribution to the free energy is 
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and Fermi functions g(ui) = (e' 3 " — 1) 1 , f% — 
( e /3(A 2 /2m-^ (T ) _ i)-i. Replacing j with x (and using 
the appropriate scattering length) gives the equivalent 
quantity in the excited state. As in the standard deriva- 
tion of the Gibbs-Duhem relation, dimensional analy- 
sis requires that the free energy density is of the form 
(30/V — A _3 /(a/A,/3/i|,/?/i|). Since in the limit of low 
down-spin density this will be proportional to the number 
of down-spins, we have that the interaction energy will 
be of the form A/fej. = E- m t = O — Oq— TS — /tfiVf — /if 
= (—3/2)80, — (l/2)ad50/da. As one tunes from the di- 
lute limit, na 3 <C 1, to unitarity, the interaction energy 
goes from E int = -250, to E int = -3/250 Formally 
one can write Eq. ^ in the limit /3/t — » — oo, 
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Where 6, defined in Eq. ((6]), is evaluated with f ^ = 0, 
implying it is only a function of v = u) + /if + /if — 
q 2 /4m, f3,q and /if. At low temperature, /3/tf ^> 1, this 
integral is poorly behaved and is best replaced by a low 
temperature expansion whose leading term is 
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where the inverse temperature p enters in the Bose 



where 6 can be calculated analytically, and is evaluated 
at v = q 2 /4m. At unitarity we find e = — 0.69/if, a rea- 
sonable approximation to the results of quantum Monte- 
Carlo [29(. Given that the lowest temperature highly 
polarized data from 0] has /3/tf = 1.25, finite tempera- 
ture effects are crucial and the asymptotic result cannot 
be compared to experiment. 

As previously anticipated, at high temperatures 
Eq. reduces to the polarization-imbalanced virial 
expansion result e = ^/2k B T\ T n^ (-36 2 /4 - ydb 2 /dy) 
(cf. jig), where the second virial coefficient is 62 = 
(l/7r)/ °°dx/(l + x 2 )e- x2 / y , with y = mk B Ta 2 /h 2 . At 
low densities (na 3 -C 1) we can analytically integrate © 
to recover e = (Airh 2 /m)ani , regardless of temperature. 

The spectra in figure [2] capture many features of the 
experiments and of more numerically demanding the- 
ories @, H[. The most successful element is that since 
we include final state interactions we are able to repro- 
duce the energy scales observed in experiment. Addi- 
tionally, as already emphasized, our intermediate tem- 
perature spectra are bimodal - however the bimodal- 
ity is much much less pronounced in our theory than 
in the experiment. By comparing with the calculation 
of Massignan, Bruun, and Stoof [8J, we can understand 
this discrepency as an artifact of neglecting the width 
of the homogeneous spectrum - due in part to the fact 
that at higher densities the initial | j.) state has overlap 
with many \x) states. The intermediate detuning spectral 
weight (for example around 10 kHz in figure [5]) will be 
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reduced by this broadening, more clearly separating the 
spectrum into two peaks. In experiments, finite probe 
duration further broadens the spectrum, and the finite 
probe size and trap anharmonicities distort the low de- 
tuning peak. 

Having understood the case n^/n^ <C 1, we now ad- 
dress the balanced situation nj_ = nj. It should be clear 
that all of our arguments about the scaling of the energy 
with density carry over to the unpolarizcd gas. Thus, 
to the extent that the homogeneous spectrum consists 
of a single peak, one expects that the trapped spectra 
should have the same qualitative features that we saw 
in figure [21 Indeed, there is very little qualitative differ- 
ence between the experimental spectra of polarized and 
unpolarized atoms in a trap [H . 

On the other hand, our quantitative calculations do not 
generalize easily to the unpolarized case. The physics of 
the unpolarized gas is more complicated than that of the 
highly polarized gas: with increasing the down-spin 
atoms occupy a larger range of momentum states, the up- 
spin Fermi sea becomes increasingly perturbed, and the 
normal state eventually becomes unstable to forming a 
superfluid. While RF spectroscopy is undoubtedly sensi- 
tive to these many-body effects, their manifestation may 
be subtle, sensitive to final state interactions, and easily 
obscured by the trap inhomogeneities. There does not 
appear to be a model independent way to use inhomoge- 
ncously broadened spectra to demonstrate superfluidity 
or extract a pairing gap. 

Finally, we mention that with minor changes our qual- 
itative arguments also apply in the simpler case of a 
weakly interacting Bose gas, such as spin polarized hy- 
drogen [3(| ■ In these gases the interaction energy is pro- 
portional to the density. At low temperature, nX^ > 1, 
the density is proportional to the chemical potential, 
while at high temperature, nXj, <C 1, the density is ex- 
ponentially small in f3[i. Following the arguments pre- 
sented here, one will therefore see a bimodal spectrum in 
a trapped Bose gas at intermediate temperatures. Such 
bimodality is indeed found [31] ■ 
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